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We explore the use of SKA to deduce the physical parameters of kiloparsec-scale jet flows in 
radio galaxies. Jets in Active Galactic Nuclei are relativistic where they are first formed, but their 
speeds and compositions change as they propagate. It has long been known that kiloparsec-scale 
jets in radio galaxies can be divided into two flavours; strong (found in powerful sources, narrow 
and terminating in compact hot-spots) and weak (found in low-luminosity sources, geometrically 
flaring, unable to form hot-spots and terminating in diffuse lobes or tails). We have developed 
methods to model AGN jets as intrinsically symmetrical, relativistic flows by fitting to deep, 
well-resolved radio images in Stokes /, Q and U. This has yielded a wealth of information about 
the brightest few weak-flavour jets. Our first key objective is to observe large samples of weak 
and transition jets at 0.1 - 0.5 arcsec resolution with SKAl-MID. This would allow us to see how 
jet propagation depends on power and environment and to quantify the energy and momentum 
input into the IGM. We will require typical noise levels of 1 /iJy/beam, and may be able to exploit 
survey imaging in some cases. Our second, more challenging, application is to determine the 
velocity fields in strong-flavour jets. Do they have very fast (T Ri 5 — 10) spines? Is there evidence 
for magnetic confinement by a toroidal field? What are their energy fluxes? This is a major 
imaging challenge for SKA2: we need resolution better than 0.05 arcsec, ideally in the 1 - 10 GHz 
frequency range, with rms noise levels 10 nJy/beam and extremely high dynamic range, imaging 
fidelity and polarization purity. 
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1. Introduction 

Jets from active galactic nuclei (AGN) are important in many areas of astrophysics: they 
extract energy from supermassive black holes, produce the most energetic photons (and perhaps 
cosmic rays) we observe, act as conveyors of ultrarelativistic particles and magnetic fields from the 
parsec-scale environments of AGN to the multi-kiloparsec scales of extended radio galaxies and 
quasars, and supply copious amounts of energy to their surroundings, thereby preventing cooling 
and profoundly affecting the evolution of massive galaxies and clusters. 

An overview of the study of AGN jets with SKA is given by Agudo et al. (2015). This paper 
covers one aspect - the physics of jets on kpc scales - in more detail, emphasising a number of key 
topics (cf. Blandford 2008), as follows. 

1. What are the velocity fields? Why do some jefs expand rapidly, increase in brighfness and 
decelerafe while ofhers remain fasl and collimated? How does Ibis relate fo fhe origin of fhe 
Fanaroff-Riley classes (Fanaroff & Riley 1974)? 

2. How does jef composifion change wifh disfance from fhe AGN? Whaf are fhe typical mass 
fluxes and enfrainmenf rales? 

3. Whaf is fhe magnelic-field lopology on kpc scales - ordered or disordered? 

4. How is jef propagalion determined by fhe external environmenl (densily, pressure, tempera- 
lure, magnelizalion)? 

5. How are jefs confined: by external gas pressure, magnelic fields, some olher mechanism or 
nol af all? 

6. Where and by whaf process are relalivislic particles acceleraled? 

7. Whaf are fhe jef energy and momenlum fluxes? 

8. Whaf effecls do fhe jefs have on Iheir surrounding IGM/ICM? Do Ihey lead lo shocks and 
heating? Whaf effecl do Ihey have on exlernal magnelic fields and hence on Ihermal Irans- 
porl? 

If has long been known lhal AGN radio jefs divide info slrong and weak ‘flavours’ (Bridle 1984). 
Slrong-flavour jefs are usually one-sided (>4:1) and well-collimaled over Iheir whole lenglhs and 
often terminate in compact hot spots. Weak-flavour jets are initially asymmetric, but tend to sym¬ 
metrize far from the nucleus. They also ‘flare’ wifh increasing disfance from fhe AGN in Iwo 
dislincl senses: geometrically (a significanl increase in apparenl spreading rale) and in brightness 
(an increase in apparenl inlensily, often following an inilial ‘gap’, or exlended region in which 
fhe radio emission is weak or undeleclable). The qualilalive piclure is lhal bolh jef flavours are 
initially relalivislic, bul weak-flavour jefs decelerafe lo subrelalivislic speeds on kpc scales, while 
slrong-flavour jefs remain highly relalivislic unlil Ihey reach fhe ouler parls of fhe sources. 

Thus far, we have derived quanlifafive models only for a small number of weak-flavour jets 
using deep VLA observalions. Strong-flavour jets are harder lo sludy as Ihey are narrower lhan 
weak-flavour jefs (Ihus harder lo resolve) and Iheir counler-jefs are fainl and difficull lo isolate 
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from surrounding lobe emission. After outlining the modelling methods and associated imaging 
requirements (Section we therefore consider in turn the use of SKAl-MID to observe large 
samples of weak-llavour jets (Section |3|) and SKA2 to model strong-flavour jets for the first time 
with our methods (Section^. 

2. Measuring the flow parameters of jets 

Although numerical (GRMHD) simulations of jet formation have made considerable progress 
recently (e.g. McKinney et al. 2013), the problem of simulating the propagation of a very light, 
relativistic, magnetized jet in three dimensions is computationally prohibitive, with poorly known 
initial conditions: no simulation can yet hope to follow a jet all the way from its formation on scales 
comparable with the gravitational radius of the central black hole to the kiloparsec scales which 
SKA can best resolve. A more empirical approach to estimation of jet flow parameters is therefore 
needed. 

One powerful method is based on the realization that kpc-scale jets are both signihcantly 
relativistic and intrinsically symmetrical, in the specihc sense that apparent asymmetries due to 
relativistic aberration are much larger than any intrinsic asymmetries. Jets emit (often highly) 
polarized synchrotron radiation. With the assumptions that the flows are (on average) axisym- 
metric, intrinsically symmetrical and stationary, the observed brightness distributions in Stokes 
I, Q and U can be modelled to ht for geometry, velocity held, intrinsic emissivity variation and 
three-dimensional magnetic held conhguration. The method was hrst described in Laing & Bridle 
(2002a); a comprehensive overview is given by Laing & Bridle (2014). The key to the method is to 
determine the velocity and inclination angle independently by comparing emission from the main 
and counter-jets in both total intensity and linear polarization. The observed polarizations in the 
approaching and receding jets are not the same because they are effectively observed at different 
angles to the line of sight in the rest frame, 6'. For antiparallel jets at angle 6 to the line of sight in 
the observed frame. 


sin 6' = [r( 1 — j8 cos 0)] ' sin 6 
sin0^j = [r(l -|-j8cos0)]^'sinS 


where subscripts j and cj refer to the approaching and receding jets, jS is the how velocity in units 
of c and F is the bulk Lorentz factor. jS and 6 cannot be determined independently from images of 
total intensity alone: the use of linear polarization is crucial. 

This method has proved to be very successful in modelling the jet hows in bright, nearby (FR I) 
radio galaxies, as we describe in Section Given the derived velocity helds (and the assumption 
of pressure equilibrium with the surroundings on large scales), we can also use the laws of con¬ 
servation of mass, momentum and energy to estimate the energy and mass huxes and the run of 
entrainment rate along the jets (Laing & Bridle 2002b). The magnetic-held topology is not com¬ 
pletely determined, however. The models show that longitudinal and toroidal components dominate 
(at small and large distances, respectively). The observed symmetry of the transverse prohles of I 
and Q tells us that the held structure is not a simple ‘grand design’ helix on these scales (Laing et 
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al. 2003), but leaves open the question of whether the toroidal component is vector-ordered. If it 
is, there should be systematic gradients in Faraday rotation across the jets (we know that thermal 
electrons must be entrained even if they are not present near the AGN). Such searches have so far 
been frustrated by large fluctuations in foreground Faraday rotation. 

The observational requirements for this technique are images of both jets with good fidelity, 
at least 5 and preferably >10 beams across their widths and noise levels low enough to detect 10% 
polarization at 4a everywhere. It must also be possible to correct the effects of Faraday rotation 
(and depolarization) to high accuracy. 

These methods have so far managed to provide answers to many of the questions posed in 
Section [T], but only for a small number of carefully selected sources. We have been restricted to jets 
which are bright and easily resolved. Thus we sample a limited (intermediate) range of luminosity, 
and can only model jets with large opening angles. We know very little about more powerful, 
highly-collimated jets, intrinsically weak sources or the systematic dependences of jet physics on 
power and environment. To do better, we need much large samples of low-luminosity sources 
(high sensitivity at moderate resolution) and an instrument capable of much higher resolution and 
very high sensitivity for the powerful sources. We explore these applications in the following two 
sections. 

3. Decelerating jets: SKAl-MID 

3.1 Current results 

Our current results on decelerating, FRI jets are summarised in Laing & Bridle (2014), and 
an example is shown in Fig. All of the jets first increase in opening angle and then recollimate 
to form conical flows at a fiducial distance ro (the flaring region) in the range 2-35 kpc from 
the AGN. At ssO.lro, the jets brighten abruptly at the onset of a high-emissivity region and we 
find an outflow speed of ?«0.8c, with a uniform transverse profile. Jet deceleration first becomes 
detectable at Ri0.2ro and the outflow often becomes slower at its edges than it is on-axis. Decel¬ 
eration continues until R^O.hro, after which the outflow speed is usually constant. The dominant 
magnetic-field component is longitudinal close to the AGN and toroidal after recollimation, but 
the field evolution is initially much slower than predicted by flux-freezing. In the flaring region, 
acceleration of ultrarelativistic particles is required to counterbalance the effects of adiabatic losses 
and account for observed X-ray synchrotron emission, but the brightness evolution of the outer jets 
is consistent with adiabatic losses alone. These results are best interpreted as effects of the interac¬ 
tion between the jets and their surroundings. The initial increase in brightness occurs in a rapidly 
falling external pressure gradient in a hot, dense, kpc-scale corona around the AGN. We interpret 
the high-emissivity region as the base of a transonic ‘spine’ and suggest that a subsonic shear layer 
starts to penetrate the flow there. Most of the resulting entrainment must occur before the jets start 
to recollimate. 

3.2 SKAl-MID observations 

While the basic process of jet deceleration in twin-jet sources has become clearer, there are still 
many unanswered questions. Almost all of the sources we have modelled so far have luminosities 
close to P ^ 10^^ ^ WHz^^ at 1.4 GHz. The next set of questions includes the following. 
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Figure 1: Observations and models of the nearby radio galaxy NGC 315 (Laing & Bridle 2014). (a) Vectors 
with magnitudes proportional to the degree of polarization, p, and directions along the apparent magnetic 
field superposed on a false-colour plot of total intensity, (b) Model fit corresponding to (a), (c) The derived 
velocity field, in units of c. 

1. How does the characteristic recollimation scale ro vary with luminosity (jet power) and ex¬ 
ternal pressure? Can we develop a way of predicting a priori which jets will decelerate and 
the form of their outer structures? 

2. What are the distributions of energy and mass flux? How do these depend on black hole mass, 
accretion rate and ISM/IGM properties? Do they correlate with estimates from the size and 
pressure of jet-blown ‘cavities’ in the IGM/ICM? Can we infer the energy and momentum 
input integrated over the source population for comparison with galaxy formation models? 

3. What causes jet emission to ‘turn on’ at the flaring point, ^ O.lro? 

4. How important is stellar mass input in decelerating the weakest jets (it can be differentiated 
from boundary-layer entrainment because no slow boundary layer is formed)? 

5. Can we decide unequivocally whether the toroidal held component is vector ordered and, if 
so, does it contribute to jet collimation on kpc scales? [To do this, we will need to measure 
Faraday rotation for a sample of sources selected to have low foreground contamination.] 

6. How important are intrinsic and environmental asymmetries? This cannot be assessed for an 
individual object, but relativistic and environmental asymmetries are expected to be uncor¬ 
related, so we can separate their effects statistically for a sufficiently large sample. 

In order to answer these questions with any statistical significance, we will need to increase the 
number of modelled sources from 10 to at least a few hundred and to identify sub-samples with 
different luminosities and environments. Our long VLA observations typically reach rms noise 
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Figure 2: VLA images of (a) the powerful quasar 3C334 (Bridle et al. 1994) and (b) the radio galaxy 
3C353 (Swain et al. 1998). Note the very narrow strong-flavour jets. 3C334 is thought to be at ^ 40° to 
the line of sight and shows a pronounced jet/counter-jet asymmetry (the counter-jet is only visible at large 
distances from the AGN). In contrast, 3C 353 is close to the plane of the sky and both jets are visible, but its 
asymmetries are dominated by the environment. 


levels between 5 and 10/xJy/beam at angular resolutions between 0.25 and 2.5 arcsec (The up¬ 
graded Karl G. Jansky VLA is allowing us to go up to a factor of 5 deeper, but only with long 
multi-configuration observations in special cases). For an integral source count N{> S) 
(appropriate for these nearby sources), we would need to increase the sensitivity by a factor of 10 
from the flux densities alone (i.e. an rms noise level of 0.5 — 1.0/f Jy/beam). Detailed trade-offs 
depend on source brightness, size and environment. For example: the need to measure low Fara¬ 
day rotations argues for lower frequencies, but we cannot tolerate a significant variation of rotation 
across the observing beam; there are systematic variations of jet size (and therefore required res¬ 
olution) with luminosity, and so on. Band 3 is therefore preferred for the measurement of low 
Faraday depths; Band 5 for jet modelling. These observations would also probe the magnetoionic 
environments of the radio galaxies in great detail (Laing et al. 2008; Guidetti et al. 2011). 

An initial estimate is that we would need fti300 observations (mostly pointed, but perhaps 
also obtained from surveys) with high image fidelity, resolution ftiO.25 arcsec and rms noise level 
^ 1 /iJy/beam in SKAl-MID, Band 5. For the combined SKAl-MID array (Dewdney et al. 2013, 
Table 9) and allowing for a factor of 2 increase in noise from weighting (Braun 2013), this would 
give an average of 50 min on-source time per observation, or ^ 240 hr in total. The observations 
would still be feasible in an Early Science phase, but longer integrations (by factors ^ 4) or 
correspondingly smaller sample sizes would be required. 

4. Fast, powerful jets: SKA2 

The strong-flavour jets found in powerful radio galaxies and quasars, particularly those associ¬ 
ated with FRII sources, are far harder to study, as they are narrow and (in the case of counter-jets) 
also very faint (e.g. Fig. §). 

We aim to distinguish between two hypotheses: (a) strong jets are only mildly relativistic (jS « 
0.5 — 0.7; F « 1.2 — 1.4) as inferred naively from the statistics of their integrated sidedness ratios 
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(Wardle & Aaron 1997; Mullin & Hardcastle 2009) or (b) they have a two-component velocity 
structure, with a much faster (j8 > 0.95) spine surrounded by a slower shear layer that dominates 
the emission except in very end-on cases. Case (b) has long been suspected (Bridle et al. 1994), 
but never proved: arguments in favour include the following. 

1. High pattern speeds are inferred from VLBI proper-motion measurements for pc-scale jets 
in flat-spectrum quasars (the end-on counterparts of sources like 3C 334; Lister et al. 2013). 
If these are representative of the bulk motion, then flow Lorentz factors up to F 40 are 
required. Deceleration of such powerful jets to mildly relativistic speeds on kpc scales cannot 
occur without disruption. 

2. If the X-rays observed in the extended jets associated with core-dominated sources are gener¬ 
ated by inverse Compton scattering of cosmic microwave background photons by relativistic 
electrons in the jet, then the bulk Lorentz factors must be large (F 10) and the jets must 
be close to the line of sight (Tavecchio et al. 2000; Celotti et al. 2001). Note, however, that 
there are arguments against this interpretation of the X-ray emission, at least in some sources 
(Hardcastle 2006; Meyer & Georganopoulos 2013). 

3. The jet/counter-jet ratio and observations of apparent superluminal speeds in the jet of M 87 
both indicate F 6 over at least part of its length (Asada et al. 2014). 

There are a number of transition cases between weak and strong-flavour jets which do not 
appear to decelerate significantly but which have counter-jets that can now be imaged with the 
VLA. The best example is NGC 6251 (Fig ^), for which observations are under way. 

Modelling of true strong-flavour jets such as those found in luminous quasars is a far more 
challenging prospect, even for SKA2. Counter-jets can currently only be imaged with something 
approaching adequate s/n in sources such as 3C 353 (Fig. Swain et al. 1998) which are very 
close to the plane of the sky (in which case our modelling technique does not work). We really 
need to image both jets in powerful sources which are reasonably inclined. As an example, if 
the straight part of the counter-jet in 3C 334 (Fig. §) has I « 10/iJy/beam at 5 GHz with a beam of 
0.35 arcsec (just below the current limit) and a width of 0.2 - 0.6 arcsec, we would need a resolution 
of Ri0.05 arcsec, and an rms noise level of lOnJy to be able to image and model the counter¬ 
jet in linear polarization. In turn, this requires a dynamic range of 10^ : 1 (peak:rms). Both 
requirements could be relaxed if adequate resolution could be obtained at a lower frequency. 

5. Summary 

We have outlined some of the ways in which SKA could contribute to the study of large-scale 
jets in radio galaxies. Extension of current work on low and intermediate luminosity sources to 
large samples is feasible for SKA 1-MID. The case of powerful jets is an interesting challenge for 
SKA2. 
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Figure 3: Images of NGC6251. (a) Overall structure (WSRT, 609MHz, from Mack et al., 1997). (b) 
Rotated image of the jets at 1665 MHz, 4 arcsec FWHM from archival VLA data, (c) Sidedness ratio R = 
/j//cj, derived from (b) and blanked (grey) where / < 3(7 in either jet. The jets in this source represent a 
transition between weak and strong flavours; they flare, but still maintain a large side-to-side ratio. The 
sidedness ratio image in panel (c) strongly suggests a substantial velocity gradient across the jets. 
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